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ABSTRACT 

We investigate how the X factor, the ratio of the molecular hydrogen column 
density (Nn 2 ) to velocity-integrated CO intensity (W), is determined by the physical 
properties of gas in model molecular clouds (MCs). The synthetic MCs are results of 
magneto-hydrodynamic simulations, including a treatment of chemistry. We perform 
radiative transfer calculations to determine the emergent CO intensity, using the large 
velocity gradient approximation for estimating the CO population levels. In order to 
understand why observations generally find cloud-average values of X= A^Gai^ 2 x 10 20 
cm~ 2 K _1 km -1 s, we focus on a model representing a typical Milky Way MC. Using 
globally integrated iVn 2 and W reproduces the limited range in X found in observations 
and a mean value X= XQ a \= 2.2 x 10 20 cm -2 K _1 km -1 s. However, we show that 
when considering limited velocity intervals, X can take on a much larger range of values 
due to CO line saturation. Thus, the X factor strongly depends on both the range 
in gas velocities, as well as the volume densities. The temperature variations within 
individual MCs do not strongly affect X, as dense gas contributes most to setting the 
X factor. For fixed velocity and density structure, gas with higher temperatures T 
has higher W, yielding X oc T^ 1 / 2 for T ~ 20 - 100 K. We demonstrate that the 
lincwidth-size scaling relationship does not influence the X factor — only the range in 
velocities is important. Clouds with larger linewidths a, regardless of the linewidth-size 
relationship, have a higher W, corresponding to a lower value of X , scaling roughly 
as X oc er -1 / 2 . The "mist" model, often invoked to explain a constant Xq^i consisting 
of optically thick cloudlets with well-separated velocities, does not accurately reflect 
the conditions in a turbulent molecular cloud. We propose that the observed cloud- 
average values of X^Xq & \ is simply a result of the limited range in Nn 2 , temperatures, 
and velocities found in Galactic MCs — a nearly constant value of X therefore does 
not require any linewidth-size relationship, or that MCs are virialized objects. Since 
gas properties likely differ (albeit even slightly) from cloud to cloud, masses derived 
through a standard value of the X factor should only be considered as a rough first 
estimate. For temperatures T ~ 10 — 20 K, velocity dispersions a ~ 1 — 6 km s -1 , and 
N H2 - 2 - 20 x 10 21 cm~ 2 , we find cloud-averaged values X~ 2 - 4 x 10 20 cm" 2 K" 1 
km -1 s for Solar-metallicity models. 

Key words: ISM:clouds - ISM: lines and bands - ISM: molecules - ISM: structure 
- line:profiles - stars:formation 



1 INTRODUCTION 

Carbon monoxide (CO), the second most abundant molecule 
in the interstellar medium (ISM), has now been observed for 
over thirty years to investigate the physical characteristics 



of the ISM. The lowest rotational levels are easily excited 
through collisions with molecular hydrogen (H2), by far the 
primary component of molecular gas in the ISM. Since H2 
is difficult to detect directly, and the 12 CO J=l-0 line oc- 
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curs at a frequency (115.27 GHz) that is readily observable 
from Earth, CO observations are well suited for probing the 
conditions of the molecular component of the ISM. 

Accurately measuring the masses and velocities of gas 
within molecular clouds (MCs) is of primary importance for 
understanding star formation. CO observations of Galactic 
and extra-galactic MCs have provided a wealth of informa- 
tion about these properties, allowing for detailed modeling 
of the star formation process. However, uncertainty remains 
about exactly how to convert observed CO emission into 
fundamental physical properties of the dominant molecular 
component of the ISM. 

The low rotational transition lines of CO are known 
to be optically thick, and therefore a considerable fraction 
of the emission from high density regions must be self- 
absorbed. Nevertheless, a strong correlation is found be- 
tween the CO intensity and the H2 column density Nh 2 ■ A 
number of methods are employed to measure Nn 2 - These 
include mass determinations using observations of 13 CO, 
which has lower optical depth than 12 CO and so may b e 
capable of tracing much of the MC gas (e.g. [Dickman 19781 ). 
Alternatively, if MCs are in virial equilibrium, the 12 CO 
linewidth may be used to estim ate the virial mass (e.g. 
lLarson|[l98ll ; [Solomon et al.lll987h . Independent mass mea- 
surements not involving CO include observations of 7-rays, 
which are produced when cosmic rays interact with the ISM. 
The molecular content is deduced when HI observations pro- 
vide information about the amount o f atomic material along 
the line-of-sight (jStrong et alJll988h . Dust-based observa- 
tions in the infrared may also provide indirect g as mass esti- 
mates, using appropriate dust-to-gas ratios (e.g. lDame et al.l 
l200ll ; iPineda et aiIl200Sl ; iLerov et alj|2009ft . 

The correlation between 12 CO (J=l-O0 intensity W 
and Nh 2 is expressed as 

X = ^(cm- 2 K- 1 km- 1 s). (1) 

In the Milky Way, observational analyses generally find 
this "X factor" to be nearly constant and ~fewxl0 20 
cm" 2 K" 1 km" 1 s, hereafter X&g, for both MCs (see e.g. 
ISolomon et ailll987l ; lYoung fc Scovilleill99ll . and references 
therein) and the lower density diffuse ISM (e.g. [Polk et all 
Il988l : iLiszt et aljl20ich . In regions with very high molecular 
densities, such as ultraluminous infrared galaxies (ULIRGs), 
X is found to be ~l-5 times lower than XGai- There is 
believed to be only a small range in the X factor among 
ULIRGs fsee lSolomon fc Vanden BoutH2005l . and references 
therein). As a result suitable X factors, depending on en- 
vironment, are commonly used for directly estimating the 
7Vh 2 (or gas mass) from CO observations. 

A key assumption in using an X factor to derive gas 
masses is that the CO line is an approximately linear tracer 
of the bulk of the MC gas. Since the CO line is opti- 
cally thick, however, it is not obvious why a linear relation- 
ship s hould hold. Exte nding the analysis of iDickman et all 
(1986). ISoiomon et al.l l|l987t ) advanced the "mist" model in 
which a MC is composed of optically thick cloudlets that 
have well separated velocities, such that the amount of self- 
absorption in the MC is small, and W is proportional to the 
total number of cloudlets along the line-of-sight, hereafter 

1 We will hereafter refer to 12 CO (J=l-0) simply as "CO." 



LoS. The assumptions implicit to this "mist" model have 
not, however, been tested with radiative transfer modeling 
in realistic MC models. 

Despite the correlations between W and A r H 2 de- 
scribed above, there are in fact signs that no univer- 
sal X factor is applicable to all sources. First, since the 
CO line is optically thick, observations of Galactic MCs 
show that beyond a threshold column density, W sat- 
urates so that the CO line no longer traces gas mass 
jLombardi et alj|2006l : [Pineda et al.ll200Sft . Observations of 
low metallicity sys t ems such as t he SMC suggest X ^ $> Xg a \ 
(llsrael et al.| 1 19861: llsraell Il997l : iBoselli et all Il997l , 120021 : 
iLerov et alj 120091 . l201ll ). Theoretically, low metallicity sys- 
tems are expected to have lower CO abundances, which 
lead to lower CO intensities and thereby larger X fac- 
tors jMaloney fc Blackll 19881 : IWolfire et all 19931 : (israelii 19971 : 
IShettv et al.ll201lf ) provided that the other properties of the 
MCs (e.g. mass, size) do not vary greatly with metallicity 
|Glover fc Mac Lowll201ll ). Further, using independent dust- 
based measures of Nn 2 leads to X factor estimates that 
differ from X computed through the "virialized" clo ud as- 
sumption (|Bolatto et al.ll200Sl ; lLe"rov et alj|2007l ,l2009). This 
discrepancy is taken as evidence that t here are large reser- 
voirs of molecular gas untraced b y CO (|Grenier et al.ll2005l : 
IPlanck Collaboration et al.ll201ll ). One explanation for such 
a situation is that in the outer regions of MCs, CO can- 
not f orm efficiently due to poor self-shielding (|Wolfire et al.l 
l201Gh . 

The discrepancies and range in the X factor estimates 
need to be understood in order to accurately interpret CO 
observations. One important issue related to the X fac- 
tor is the dynamical state of MCs. Two well known scal- 
ing relationships for MCs have been established in large 
part due to CO o bservations: t he mass-size and linewidth- 
size relationships l|Larsonlll98ll ). often referred to colloqui- 
ally as "Larson's Laws." Assuming constant X, the MC 
masses M are found to be related to the projected size R 
through a power law with approximate index 2, M oc R 2 . 
Observed linewidths for MCs follow a power law relation- 
ship with projected size a oc R 1 ^ 2 . Taken together, the in- 
terpretation of these relationships is that clouds are (ap- 
proxima tely) in virial equilibrium, or "viria l ized," so a 2 ~ 
GM/R (lLarsonlll98ll: IDickman et al.lll986l : ISolomon etafl 
1 19871 : iMvers fc Goodmanl 119881 ). Alternatively, if MCs are 
virialized such that Nh 2 ~ M/(fiR 2 ) ~ a 2 /(^GR), then 
X ~ a 2 /(pGRW). If a 2 /R and W vary little over the pop- 
ulation of MCs, then X would have a uniform value. 

An open question is whether the "mist" model is appli- 
cable to turbulent media, as turbulence is now considered 
a dominant factor controlling the dynamics of M Cs (e.g. 
iMac Low fc Klcsscn 2004; McK ee fc Ostrikeill200l and ref- 
erences therein). In the first publication in this series 
|Shettv et al.ll201ll . hereafter Paper I), we investigated how 
well CO can trace the underlying molecular gas using radia- 
tive transfer calculations on magnetohydrodynamic (MHD) 
models of MCs, which include a treatment of chemistry. 
We focused on velocity integrated CO intensity, and com- 
pared probability distribution functions (PDFs) of W, Nu 2 , 
and CO column density Nco- We also assessed the X fac- 
tor, again only considering velocity integrated intensities. 
We showed that even though X may vary between differ- 
ent LoSs through a given MC, the cloud-averaged intensity 
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produces XmXaaX within a factor of 3, from various mod- 
els with different metallicities Z/Zq = 0.3 — 1 and densities 
no = 100 — 300cm -3 . In this work, we focus primarily on 
the Milky Way MC model from Paper I, and investigate 
the properties that affect the derived X factor. Using radia- 
tive transfer calculations of turbulent chemo-MHD models, 
we perform a rigorous theoretical investigation of the qual- 
itative models proposed to explain the observed X factor. 
We directly modify the cloud characteristics, such as tem- 
perature, density, and velocity, and recompute the X factor 
to understand its dependence on those parameters. Two of 
our goals are to understand why the X factor is roughly 
constant for a range of systems, including Milky Way field 
GMCs, and to assess whether the "mist" model is applicable 
to turbulent MCs. 

This paper is organized as follows. In the next section, 
we provide an overview of the estimated X factor in various 
environments. We also discuss the results from Paper 1 in 
the context of observationally derived values. In Section[3]we 
review our method of modeling CO emission from turbulent 
MCs, and discuss some properties of the main Milky Way 
MC model. In Section [4] after discussing how the X factor 
is measured, we compare how X depends on the dynamic, 
chemical, and thermal structure of the model MC. We in- 
vestigate in detail the dependence of X on various cloud 
characteristics, such as temperature, H2 and CO densities, 
as well as the velocities. In Section [5] we compare our re- 
sults to previous observational and theoretical efforts, and 
offer a new explanation for observed trends. Section [6] sum- 
marizes our interpretation of the observed X~XG a iand our 
conclusions regarding the parameter dependences of X. 



2 OVERVIEW OF THE X FACTOR 

As alluded to in the previous section, the X factor is found 
to be nearly constant when considering a specific class of 
sources, such as Galactic MCs, or a different value for 
ULIRGs. Yet, analysis of the full spectrum of molecular en- 
vironments generally portray a trend of decreasing X factor 
with increasing molecular surface density. Figure Q] shows a 
compilation of observationally in ferred X factors from vari- 
ous systems l|Tacconi et al ] l2ooah . 

In Paper I, we analyzed the X factor in various models 
with different metallicities Z/7jq = 0.1 - 1 and densities no = 
100 — 1000 crn~ 3 . Figure [2] shows the relation between the 
X factor and surface density E gas or Nh 2 from the various 
models discussed in Paper l|f| The points show the mean X 
factor averaged in bins of E gas , and the large points show the 
cloud-averaged X factor. At intermediate densities £ gas ~50- 
200 Mq pc -2 , the cloud- averaged X factors for all but the 
very low Z=0.1 model are ~ 2 — 5 x 10 20 cm~ 2 K -1 km -1 
s « XgbI- 

Both Figs. [T] and [2] show that X is larger than Xaal 
in low metallicity systems. As explained in Paper I, large 
X factor values at low metallicities and densities are pri- 
marily due to the low W, the integrated CO intensity, 
relative to A r H 2 (Eqn. [TJ. Systems with lower abundances 



The conversion between Nn 2 an d £ ga s includes a factor of 1.4 
to account for the mass of helium. 



of CO will of course have lower W. The formation of 
CO is highly dependent on the MC density, metallicity, 
and strength of the backgr o und UV radiation field (see 
also iMaloney fc Black! Il988l; Ivan Dishoeck fe Black! Il988l ; 



iGlover et al.N2010l ; Idover fe Mac Lowll201ll ). However. H 2 
formation is not as sensitive to these properties, due to its 
ability to effectively self-shield. Thus, the relative abundance 
of CO compared to H2 can vary significantly within a MC, 
leading to a wide range in the X factor, even though the 
(emission-weighted) average of many different clouds may 
all result in a value ~AGai- 

Note that the low surface density regimes in Figs. [T] and 
[2] are not directly comparable. The low surface densities in 
our numerical models correspond to diffuse regions in MCs 
with size ~10 pc. In Figure [T] the low surface density cases 
correspond to observations with low beam filling factors. 
These regions can have sizes ~ kpc (or larger). Therefore, 
the objects in the two figures may differ in size. 

At large surface densities (£ gas ~ 200 Mq pc -2 ), our 
turbulent models show an increase in the X factor, whereas 
the observations suggest a decreasing X factor. The increase 
of X with increasing surface density in our models occurs 
because in this regime the CO line is saturated, so that W 
remains constant even as Ah 2 increases. 

Our models only account for the effects of MHD, ther- 
modynamics, and chemical evolution. The sources at high 
surface densities in Fig. [T] are either the Milky Way cen- 
ter or galaxies undergoing intense star formation activity 
(LIRGs and ULIRGs) . The heating associated with star for- 
mation, as well as the large-scale rotation and turbulence of 
the ISM in these galaxies, are not captured by our models. 
These processes likely contribute to setting the X factor in 
such environments, and may be responsible for the observed 
trends in Figure Q] at high E gas . 

In our current investigation of the X factor, we aim 
to understand which MC properties are responsible for 
A~AG a i in the 50 - 200 Mq pc -2 range. To carry out our 
analysis, we perform simple experiments by manually chang- 
ing a number of physical parameters of our models. We then 
recompute the X factor, and assess which of the modified 
parameters most affect the resulting value of X. 

2.1 Definition of the X factor 

We motivate our choice of the parameters to be modified 
by the definition of the X factor. When the CO intensity 
is expressed in units of the Rayleigh- Jeans "brightness tem- 
perature," Tb then 



W - 



J T B dv (Kkms -1 ) 



(2) 



The H2 column density is simply the volume density inte- 
grated over the LoS ds: 



Nh 2 = / n H2 ds 



Taking Equations I H3I together. 
J nH 2 ds 



X 



f T B dv 



(3) 



(4) 



This indicates that the X factor is explicitly dependent on 
three quantities: the column density of H2, the peak CO 
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Figu re 1. Compilat ion of estimated X factors from a range of systems, shown as a function of surface density. Figure reproduced from 
iTacconi et"aH j2008l) . 




Figure 2. Mean X factor in bins of gas surface density £ gas (bottom axis) or Nn 2 (top axis) for 5 models. The X factor is averaged in 
different £ gas bins. The value of X is plotted on the midpoint value of £ gas of each bin. Each model is identified by different colors and 
symbols (and labeled in the legend). The large symbols shows the global (emission- weighted) mean X factor and mean E gas from each 
model. 



intensity, and the range in velocities. Due to the coupling 
between hydrodynamics, thermodynamics, and chemistry, 
Tb is also dependent on the velocity and density (as well as 
the kinetic temperature). We aim to understand the relative 
contribution of each of these three properties of the MC. 
After assessing the X factor from the original Milky Way 
MC model, we alter one of these properties at a time, while 
keeping the others fixed, and recompute the X factor. In this 
manner, we can identify the most important cloud properties 
responsible for setting the X factor. 



3 MODELING METHOD 

3.1 Numerical magnetohydrodynamics, chemistry, 
and radiative transfer 

To investigate how MC characteristics affect the X fac- 
tor, we analyze magnetohydrodynamic (MHD) models of 
molecular clouds that include a time-dependent treatment 
of chemistry. We perform radiative transfer calculations on 
these numerical models, in order to solve for the CO level 
populations and compute the emergent CO intensity. The 
ratio of the H2 column density to the emergent CO inten- 
sity then gives the X factor (Eqn. [1]). 

The MHD grid-based models follow the evolution of 
an initially fully atomic medium with constant density in 
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a (20pc) 3 periodic box. Thermodynamics is coupled with 
chemistry to follow the formation and destruction of 32 
chemical species, including Eh and CO, through 218 chem- 
ical reactions. Emission from CO and C + are the primary 
cooling mechanisms in the dense and diffuse regions, respec- 
tively. Additionally, a constant background UV radiation 
field is included, which can photodissociate molecules in re- 
gions with insufficient shielding. The photoelectric effect is 
responsible for most of the heating in the diffuse regions, 
and in more dense regions heating is primarily due to cos- 
mic ray interactions. Turbulence is driven on large scales 
(with wavenumbers 1 $C k ^ 2) to maintain an approxi- 
mately constant root mean square (3D) velocity « rms . 

For a thorough de scription of the mo d eling m ethod, 
we refer the reade r to Idover fc Mac Low! lj2007al lbh and 
I Glover et alJ l)2010l ). In this work, we primarily consider the 
fiducial model chosen to match typical Milky Way MC con- 
ditions, which has an initial hydrogen nuclei density no=300 
cm -3 , a metallicity Z/Zq= 1, a background UV radiation 
field!] 2.7 x 10 erg cm s , and a time averaged turbu- 
lent velocity u r ms=5 km s _1 . The magnetic field is initially 
oriented along the z-axis, with magnitude 1.95 uG. 

We use the radiative transfer code RADMC- 3E0 to cal- 
culate the emergent CO intensity. The level populations in 
each zone are calcu lated through the Sobolev approxima- 
tion (|Sobolevlll957h . which uses velocity gradients across 
the faces of each grid zone to estimate photon escape prob- 
abilities. We emplo y the Einstein and collisional rate coeffi- 
cients estimated by Yang et al. (120101) and provided in the 
LAMDA database ( Schoier et alj l2005). A full description 
of the radiative transfer calculation is provided in Paper I. 

The radiative transfer calculations provide the CO in- 
tensities at each LoS position (x,y) at a given frequency v 
(or velocity) bin, I v (x, y). Besides the viewing geometry and 
the spectral resolution of the synthetic observation, the only 
other user defined parameter required to perform the calcu- 
lation is the microturbulent velocity « m trb (see Eqns. 2-7 in 
Paper I). For our fiducial model, we use « mtr b= 0.5 km s _1 . 
In the Appendix, we demonstrate that W or X does not 
strongly depend on this choice, for w m t r bS 0.25 — 0.75 km 

Since we are interested in emission from all the gas in 
the simulation volume, the spectral channels span the full 
range in gas velocities. The position-position- velocity (PPV) 
cube provided by the radiative transfer calculations has a 
spatial resolution of 0.08 pc or 0.16 pc, corresponding to the 
extent of the zones in the 256 3 or 128 3 simulations, respec- 
tively, and a spectral resolution of 0.06 km s _1 . 



3.2 The fiducial Milky Way GMC model 

Figure [3] shows JVh 2 and W of the fiducial Milky Way model 
MC, n300. Though the overall morphology of Nh 2 is ev- 
ident in the CO image, there are some stark differences. 
Most notably, the brightest regions in the CO map (near 
the bottom) do not correspond to the region with the high- 
est column density (near the top right). This discrepancy 



Table 1. Characteristics of Standard "Milky Way GMC" Model 
(n300) 



Property 



Value 



Box Size L 
Initial Atomic Density no 

{ n H 2 )vol 

(T) 
(T) vo l 

&v,los,mass — [( v f os ) mass] 2 
a v,los,vol = [( v f s)vol} 2 

{X) re f 



20.0 pc 
300.0 cm- 3 

1098.0 cm" 3 
145.9 cm -3 

19.8 K 
51.7 K 

2.4 km s" 1 
2.4 km s" 1 



1.9 xlfr 



cm 



' 2 K" 1 km" 1 



3 Thi s value i s 1.7G p as determined bv lDraind il978f) . where Go 
is the lHabind jl96Sf > field. 

4 www.ita.uni-heidelberg.de/~dmlemond/software/radmc-3d/ 



arises due to the high optical depth in the CO line. The dif- 
ferences between the observed emission and the underlying 
gaseous properties of these and other models are described 
in more detail in Paper I. 

Table [T] lists the mass and volume weighted H2 
number densities, temperatures, and LoS rms velocitiefl 
{ v los) X ^ 2 = a v,ios, along with the initial density no and 
box size L of model n300. Mass and volume weighted 
quantities are defined as (/} maS s = j ' pdV/ pdV and 
if) vol = Y,f dv /J2 dV ' respectivelyljThe last row of Ta- 
ble [1] shows a reference value of the X factor: {X) re f = 
(n H2 )voiL/[(T) 

mass&v ,los \ ■ 

As discussed in Paper I and Idover et all l|2010f ), the 
internal properties of the fiducial Milky Way cloud model 
take on a range on values. For example the CO abundance 
relative to H2, fco, can vary from ~ 10 -9 to 10 -4 . Sim- 
ilarly, the temperature, density, and velocity also take on 
a wide range of values. Throughout our analysis, besides 
considering CO emission from the original model, we also 
consider models for which some relevant physical character- 
istics are modified, such as those with constant temperature 
or CO abundances. In this manner, we can quantitatively 
assess the effect of the various physical characteristics on 
the emergent CO intensity, and ultimately on the X factor. 



4 RESULTS 

4.1 Measuring the X factor 

We begin by calculating the X factor through its traditional 
definition given by Equation [T] Figure f3] shows the mean X 
factor in bins of Nn 2 as a function of Nh 2 from all LoSs 
through the fiducial modeQ. As explained in Paper I, at the 
highest densities W does not increase with increasing Nn 2 

5 Since the viewing angle is along the z-axis, vi os = v z . Due to 
large-scale stochastic variations, the dispersion in v z is not exactly 
equal to 

f rms 

6 ( n H 2 )vol is not exactly 150 cm~ 3 because a sma ll fraction of 
hydrogen remains atomic (Glover & Mac Low 201lh . 

7 Though the figures presented here only show the results from 
one viewing angle (along the £-axis) we have verified that the 
results are not sensitive to any chosen viewing direction. 
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Figure 3. a) Column density Nn 2 and b) integrated CO intensity of the Milky Way model MC. 



due to the saturation of the CO line, resulting in XocNn 2 - 
Nevertheless, the mean X factor only varies between (1.5 - 
4)xl0 20 cm -2 KT 1 km -1 s for logiV H2 = 21-22.5. Given this 
limited range and the extent of the error bars, a constant 
value at its emission-weighted mean valu^fl (shown by the 
large filled circle in Fig. [4]) adequately describes the X fac- 
tor for this model. This mean X factor value {X)=2.2x 10 20 
is in good agreement with the reference value provided in Ta- 
ble [I] (X) = 1.9xl0 20 . We now consider the terms 7Vh 2 and 
W in detail, and their relationship to the physical properties 
of the model. 

Instead of integrating Equation [3] over the whole cloud, 
one could only use a limited range in velocity (or frequency, 
u). 



H 2 ,w 



Uh 2 dv. 



(5) 



If JVh 2 .k is computed along the same LoS and in velocity bins 
that match the spectral channels of the CO observation, then 
a column density cube of Nb. 2 ,v can be constructed which 
has the same configuration as the observed CO position- 
position-velocity (PPV) cube. 

We can now define the X factor at each PPV location, 
X v , associated with column density N# 3>u , intensity Tb and 
adopted channel width Av (in our case, 0.06 km s" 1 ): 



X v 



Nh 2 



(6) 



T B Av' 

Figure [5] shows how this "3D X factor" depends on the 
column density for the Milky Way model cloud. Here, Nh 2 ,v 
extends to much lower values than the values of Nn 2 shown 



in Figure 13] since the volume densities are integrated over a 
more limited region (in velocity space). 

The emission weighted mean "3D X factor" is equal 
to X factor computed in its usual way, shown in Figure U 
However, at high densities, X v is significantly larger that 
X, reaching values XL 10 22 cm -2 K _1 km -1 s, whereas the 
2D X factor is always < 10 21 cm" 2 K" 1 km" 1 s. From 
Figure [5] it is clear that the CO line is saturated at column 
densities (per frequency bin) above the mean column density 
Nn 2 , v = 3.5 x 10 19 cm -2 , as indicated by the iV H2 ,„/ (TbAv) 
line, where Tb is the mean brightness temperature in PPV 
regions where Nn 2 ,v> Nn 2 ,v 

The discrepancy at high column densities between X 
and X v is due entirely to the differences between W and 
TbAv. On a given LoS, W is the summation of TbAv 
through all velocities in the PPV cube. Figure [5] shows 
these quantities as a function of Nh 2 or Nn 2 ,v, respectively. 
Clearly, W continues to increase beyond Nu 2 ^ 10 21 cm -2 , 
whereas TbAv is saturated. Since Av is constant, this means 
that Tb saturates at a value ~ 13 K0 

Line saturation occurs when the optical depth at a given 
frequency v, r v £ 1. The optical depth along a LoS of length 
s is given by 



(7) 



number density in level i, and the Einstein (stimulated ab- 
sorption/emission) coefficients. The normalized line profile 
function 0„ depends on the LoS and microturbulent veloci- 



8 We denote any weighted-mean value with {...), whereas simple 
averages with an "overline" (e.g. W .) 



9 This saturation brightness temperature is estimated from Fig- 
ure [6] which shows the mean values of TbAv in bins of Nn 2 „. 
The maximum Tb can reach values ^ 30 K. 
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Figure 4. X factor from the Milky Way model molecular cloud. 
The open circles and error bars indicate the mean and stan- 
dard deviation of the X factor in bins of N^ 2 . The solid cir- 
cle shows emission weighted mean X factor of the whole model 
((X)=2.2x 10 20 cm~ 2 K _1 km" 1 s) at the mean column den- 
sity (Nh 2 = 9-0 X 10 21 cm -2 , corresponding to S gas = 202 M 
PC" 2 ). 

ties, as well as the kinetic temperature (see Section 2.2 of Pa- 
per I for more details). Together, the density, temperature, 
and velocity structure of the molecular cloud determines the 
optical depth, and therefore where the line is saturated. 

The mean and minimum value of r„, in bins of Nh 2 ,v, 
are also shown in Figure [6] with the scale given on the right 
ordinate. Though TbAv does not vary much in the range 
of the high column densities shown in Figure [6l a trend 
of increasing TbAv with column density is apparent until 
the Nh 2 , v ~ 10 20 - 5 cm -2 . At column densities Nh 2 ,v~ 10 21 
cm -2 , all values of r„> 1, indicating complete saturation of 
T B Av. 

That the velocity-integrated brightness temperature W 
still increases beyond Nu 2 ~ 10 20 ' 5 cm -2 , where the line is 
saturated, indicates that numerous optically thick regions at 
different v contribute to W. This results in the saturation of 
the integrated intensities occurring at higher column densi- 
ties than the saturation of the intensities in the PPV cube, 
as evident at the highest -/Vh 2 in Figure [(J] 

One consequence of the line saturation described above 
is that the amount of gas that is untraced depends not only 
on the density, but also on the velocity. For instance, a 
"cloudlet" or parcel of dense gas with similar velocity as 
another parcel along a LoS, but separated in space, would 
not contribute extra flux to Tb, or the W map. On the 
other hand, a "cloudlet" with much different velocity along 
the LoS would be traced, detected as higher Tb at a dif- 
ferent location in the velocity axis of the spectrum, and 
thus contribute to the integrated intensity W . Consequently, 
LoSs with (dense) gas spanning a larger range in velocities 



Figure 5. X v vs. iVH 2 ,i/ from the Milky Way model molecu- 
lar cloud. The mean and standard deviation of X v in bins of 
iVjj 2i „ are shown as circles with error bars. The solid square shows 
emission weighted mean X v of the whole model ((X„)=2.2x 10 2t) 
cm -2 K — 1 km -1 s) at the mean column density (Nn 2tl/ = 
3.5 x 10 19 cm" 2 ). Line shows N H2:V / (Tb~Av) at Nn 2}l/ > Nn 2 , v . 

will have larger Ws. This concept of optically thick dense 
"cloudlets" is t he ba sis of the "mist" model put forth by 
ISolomon et alj l|l987l ) to explain the uniformity in the X 
factor in Galactic GMCs. We return to the applicability of 
the "mist" model in Section 4.5. 

The differences in the relationships between W and 
TbAv with column density (Fig. [6} lead to the variations 
between X and X v (Figs. H and ©. The (2D) X and (3D) 
X v distributions are shown in Figure [7] Though the peak 
value of the two distributions are similar, there is a much 
larger distribution in X v . In considering synthetic obser- 
vations with different spectral resolution, we find that the 
range in X v depends on Av, but that (X„)=2.2x 10 20 cm -2 
K" 1 km -1 s consistently. As discussed, there are lower col- 
umn densities Nn 2 ,u in a PPV cube compared to Nn 2 . For 
such densities, the line is often optically thin, and there is 
a range of possible emergent intensities for a given column 
density (Fig. [BJ. This results in a larger range in X v com- 
pared to X. 

Taken together, the limited range in the 2D X factor 
shown in Figure [4] and [7] only occurs when the CO inten- 
sities are integrated over all velocities, and densities over 
the whole LoS. Effectively, a limited range in the X factor 
only results when combining the detailed velocity, as well as 
density, structure along the LoS. This is indicative that the 
velocity range plays a crucial role in setting the X factor. 

Of course, integrating the observed brightness temper- 
atures over a limited range in velocity is not practical, as 
obtaining Nn 2 ,v is not readily feasible in observational data 
sets. Nevertheless, such an analysis has underscored that the 
X factor depends on the total velocity width of the cloud, as 
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Figure 6. Intensity (left ordinate) and optical depth (right ordi- 
nate) plotted against column densities. Triangles show the mean 
integrated CO intensity W as a function of total column den- 
sity /Vh 2 (2D). Stars show the analogous relationship from (3D) 
PPV cubes, TgAv vs Nn 2V . Filled circles show the mean optical 
depth t„, and the open circles show the minimum value of t„ in 
the Nft 2tL/ bins. Dashed line corresponds to t v =X. 



X = N H2 /W 




20 22 
og (X) (crrT 2 K~' km" 1 s) 



Figure 7. PDFs of X factors shown in Figures [4] and [5] 



explicitly evident in Equation [4] But, how does the detailed 
velocity structure, such as its relationship with the size of 
emitting regions, affect the X factor? Further, what is the 
relative contribution of J dv compared to the other terms 
uh 2 an d Tb in Equation [4p To address these questions, we 
investigate in detail the role of each quantity in determin- 
ing the nature of the CO emission, and their relationship to 
intrinsic cloud properties. 



4.2 X factor dependence on temperature 

Due to the thermodynamics in the chemical-MHD model, 
the gas has a range of temperatures. As indicated in Ta- 
ble [T] the model MC has a (volume-weighted) mean T ~ 
50 K, with a dispersion a ~ 44 K. This range in tempera- 
tures will also result in a range in the X factor, since the 
observed brightness temperature depends on the gas kinetic 
temperature (Eqn. [2]| . 

To test the sensitivity of the X factor to the tempera- 
ture distribution, we artificially set all temperatures in the 
model to a constant value, and then perform the radiative 
transfer calculations. The resulting maps are then compared 
with the original CO map. Any discrepancies can be at- 
tributed to the differences in temperatures. 

This is illustrated in Figure [5] where (a) shows the orig- 
inal distribution of X, and (b) - (d) show the X factors for 
constant temperatures of 25, 50, and 100 K. The histograms 
from the model with the original temperatures are provided 
in each panel as dashed lines. 

The constant temperature T = 20 K model, equal to 
the mass- weighted temperature (see Table (TJ, provides an X 
factor distribution that is very similar to the original model 
(Figs. [Hp). This is expected since most of the CO is located 
in dense regions where T ^ 20 K. The relationship between 
X, Nn 2 , and iVco is also rather comparable to the original 
relationship. Note that for both the original and constant 
T=20 K models, the highest and lowest values of Nh 2 have 
large X. (cf. Fig. [4}. The differences between the 20 K and 
50 K model suggests that the X factor is more sensitive to 
the mass- weighted rather than the volume- weighted temper- 
atures. 

Figure [5] shows the mean X in bins of Nn 2 from the 
original model (as in Fig. |3} , along with two reference val- 
ues obtained directly from the simulation. The red stars 
show X computed using the mass- weighted LoS velocity and 
global mass-weighted temperature. The blue triangles show 
the corresponding X from the the volume-weighted quanti- 
ties. Along LoSs with column densities lower than 10 21 ' 8 
cm -2 , the CO fraction decreases, thereby decreasing W (see 
Paper I and Section I4.3.2J1 . Thus, the simple relationships 
between X, T, and vi os cannot reproduce the X factor at low 
column densities. At column densities Si 10 21 ' 8 cm -2 , how- 
ever, the X factor computed using the mass-weighted tem- 
perature agrees fairly well with the original X factor. The 
X factor computed from the volume-weighted temperature 
systematically underestimates the X factor. These trends 
indicate that the densest regions contribute most to setting 
W, and thereby the X factor. 

As indicated in Figure [8] clouds with higher tempera- 
tures produce lower X factors, as would expected from Eqn. 
[TI21 From (X) = 1.8 x 10 20 , 1.1 x 10 20 , and 7.9 x xl0 19 cm" 2 
K -1 km" 1 s at T = 20, 50, and 100 K, respectively, we find 
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Figure 8. The variation of the X factor with gas temperature. The position of each point shows the relationship between Nn 2 and Nqq. 
The color of each point identifies the value of the X factor, computed through Eqn. [l] with the color scale and distribution given in the 
inset plots. The emission-weighted mean value (X) is indicated in each plot, a) Original model n300, b) model n300, but with constant 
T=20 K, c) constant T=50 K, and d) constant T=100. In b), c), and d), X factor distribution from original simulation in a) is shown 
as thin dashed histogram. 
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Figure 9. X factor (black circles) from the Milky Way model 
molecular cloud, as in Fig. [J] The emission- weighted mean X at 
the mean column density is shown by the large solid circle. Red 
stars show the the X factor computed by Nn 2 / (T) maas cr v i os , 
and blue triangles show Nu 2 /{T) vo ia V: i os . <T) m(lss and (T) vol 
are the global mass- weighted and volume- weighted temperatures, 
given in Table[T] <r„ ios is the mass-weighted rms velocity (which 
is equal to the volume- weighted rms velocity) along each LoS. 



an approximate scaling {X) oc T~ ' 5 . The range in the X 
factors in the modified-temperature models is, however, sim- 
ilar to the original model. This indicates that the range in 
temperatures in model n300 is not responsible for the range 
in the X factor. As a result, it is either the range in densi- 
ties or the velocities which produce the distribution in the 
X factor. 

We note that an increase in temperature to 100 K 
is likely accompanied by other variations in the molecular 
cloud structure, such as a decrease in riH 2 an d possibly fco 
if the high temperature is due to a high star formation rate. 
As such, models for which the temperatures are manually 
scaled to high values are likely omitting real environmental 
effects that would affect the X - T relationship. 



4.3 X factor dependence on density 

We turn our attention to the X factor dependence on den- 
sity. Any given cloud will have a range in volume and column 
densities, in both Eh and CO. If W depends linearly on Ah 2 , 
then there would certainly be a constant X factor along all 
LoSs. However, as discussed in Paper I, W generally does 
not have a straightforward scaling with the column density 
Aco, especially at large Aco due to the high opacity of CO 
(see Section |4Tj. Further, Aco does not directly trace Ah 2 . 
This lack of correlation between W and Ah 2 results in a 
distribution of X along different LoSs for a given MC. We 
focus on how the volume density and CO abundance affect 
the X factor. 
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Figure 10. X factor vs. Nn 2 from model nlOOO, which has box 
size (20 pc) 3 . The line shows X = N H2 /W, with W = 70 K km 
s _1 . The large circle shows (X)=3.8xl0 20 cm -2 K _1 km -1 s at 
N h 2 = 3.0 x 10 22 cm" 2 , or S gas = 672 M pc" 2 . 

4-3.1 X factor dependence on volume density 

To test whether the X factor is more sensitive to tih 2 or to 
Nn 2 , we consider a model with a higher volume density than 
model n300. This model (nlOOO), only differs from model 
n300 in its initial atomic density of 1000 cm -3 . As described 
in Paper I, most of the carbon in this model is incorporated 
into CO, and so there is a very good correlation between 
Nco and Nn 2 (see Fig. 5c in Paper I). 

As explained in Paper I, for this high density model 
the effect of line saturation is clearly apparent in the X 
factor computed through Equation [1] Figure [10] shows the 
relationship between X and Nh 2 for model nlOOO. The line 
is not a fit, but rather the relationship X = Nh 2 /W, where 
W — 67 K km s _1 is the mean integrated intensity from 
this model. We see that X invariably increases with JVh 2 > as 
would be expected if the line were saturated. 

The X-Ah 2 relationship for model n300 (Fig. [4j is 
quite different from the relationship from model nlOOO, due 
partly to the larger maximum column density in model 
nlOOO and lower minimum column density in model n300. 
We investigate whether the total column density, which oc- 
curs explicitly in Equation [T] is responsible for the increase 
in X at large A r H 2 in model nlOOO. We perform the radiative 
transfer calculations on a model which is similar to nlOOO, 
but whose box length is decreased by a factor of ~3 to 6 pc. 
In this nl000-T6 model, the lower column densities allow 
the background UV radiation to penetrate further into the 
cloud, leading to more CO photodissociation. Nevertheless, 
this model should have similar H2 column densities to the 
n300 model, since H2 can self-shield more efficiently. On the 
other hand, the volume densities in nl000-L6 will be very 
similar to those in the nlOOO model. 




21.0 21.5 22.0 22.5 
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Figure 11. X factor vs. Nn 2 from model nl000-L6, which has 
initial density 1000 cm -3 and box size (6 pc) 3 . Line shows 
X = N H2 /W, with W = 49 K km s _1 . The large circle shows 
(X)=1.9xl0 20 cm" 2 K- 1 km- 1 s at ~N H2 = 9.1 X 10 21 cm" 2 
corresponding to S gas = 204 Mq pc~ 2 . 

Figure [TT1 shows the X— Nh 2 relation for this nl000-T6 
model. The line shows X = Nh 2 /W, with W = 49 K km s _1 
for this model. Here, (X) is similar to the n300 model. The 
X— Nh 2 relation is not as well reproduced by the Nh 2 /W 
line as model nlOOO in Figure [TOl Nevertheless, the X factor 
clearly increases with iVH 2 • Compared to model n300 (Fig. 
gj), X is similar for 7Vh 2 ~ 10 22 - 10 22 ' 5 cm~ 2 , but there 
are significant differences at lower -/Vh 2 , especially ^ 10 21 ' 5 . 
Thus, the X factor, even when calculated using the total 
column density (Eqn.\Tj), depends on the volume density nu 2 , 
rather than just the column density Nh 2 ■ Note, however, that 
when averaged over all lines of sight, (X) decreases by only 
15% even though no increases by a factor 3.3. 

4-3.2 X factor dependence on CO abundance 

Besides containing a range in densities, there is also a range 
in CO abundances at a given H2 density, as is evident in 
the co lumn density relat ionships in Fig. [S] (see also Paper 
I, and ldover et al.ll2010h . In order to investigate how this 
distribution contributes to the X factor, we reset the CO 
abundance with a constant /co= "co/"h 2 ratio, and then 
perform the radiative transfer calculations. 

Figure [12] shows the X factor - Ah 2 relation for (a) 
the original model, as well as those with (b) constant fco = 
10 -4 , (c) 10 -5 , and (d) 10~ 6 . Also shown in each panel is the 
line corresponding to the X = Nh 2 /W relation, indicating 
complete saturation of the CO line. 

With the original abundances, the X = Nh 2 /W re- 
lationship is close to the model values only at the highest 
column densities. At densities ^, 3xl0 21 cm -2 , the over- 
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Figure 12. The X factor— iVji 2 relationship from model n300, as in Figure[4](a), but with the CO number density reset to be (b) 10 4 , 
(c) 10~ 5 , and (d) 10" 6 X njj2- Line in each plot shows X = Nx 2 /W relation. 



plotted line underestimates X. In this regime, CO emission 
is not saturated. Figure [T2b is analogous to the limiting sce- 
nario where all the atomic carbon and oxygen is converted 
to CO, so that /co = 10" 4 . For this model, the X = N H2 /W 
relation is similar to the data above -/Vh 2 ~ 10 21 ' 5 cm" 2 , 
indicating that the the CO line is nearly fully saturated. In 
Figure [T^t -d. with lower fco, line saturation becomes less 
and less important, especially at lower column densities, fi- 
nally resulting in a constant X and Nu 2 relation for fco = 
10" 6 . 

In general, the CO line becomes saturated in regions 
with the highest CO abundance. With constant fco~ 10~ 4 , 
X increases with increasing Nn 2 everywhere, and the CO 
line is completely saturated at column densities ^ 10 21 ' 5 
cm" 2 (Fig. \l2h), resulting Xoc Nu 2 . At lower CO abun- 
dances, W increases with increasing column density. This 
results in a shallower slope in the X - Nu 2 relation. At very 
low fco= 10~ 6 , W is directly proportional to Nu 2 , so that 
X is constant at all Nh 2 ■ Taken together, clouds with both 
low and high CO abundances will tend to have a more lim- 
ited range in X than a cloud with only large /co, which 
would have Xoc Nn 2 ■ 

Notice that in Figure fl2ti . X= 10 21 cm" 2 K" 1 km" 1 
s, for the model with CO abundance /co = 10~ 6 . This is 
the quoted abundance in dif fuse Milky Way g as observed 
by iBurgh et all d2007T ) and iLiszt et all (|2010T ). However, 
iLiszt et all (|2010l 'l find X« 3 x 10 20 » X Ga i, which is a 
factor ~5 lower than the resulting value in Fig. I12H . The 
discrepancy is likely due to the combination of the differ- 
ences in temperature and linewidths between the "diffuse" 
MCs and more massive giant MCs. The diffuse ISM has a 
higher temperature, up to ~100 K. As we demonstrate in 
Section [4.21 such high temperatures may account for a sig- 
nificant fraction of the difference. Further, observations of 
low column density LoSs probably trace a larger volume of 



the Galaxy, thereby including gas with a wider range in ve- 
locities than those found in MCs. As we discuss in the next 
section, larger velocities may lead to lower X factor values. 
Nevertheless, note that the discrepancy may be partly due 
to the fact that in this numerical experiment, we artificially 
fix the CO abundance. The metallicity and self-shielding are 
not self-consistently tracked in these experiments. 



4.4 X factor dependence on velocity 

In order to test the sensitivity of the X factor to the velocity 
vios or the integral over dv appearing in Equation [4] we con- 
sider MC models with different velocity fields. In the MHD 
simulations, turbulence is generated by continuously driving 
the gas velocities with uniform power between wavenumbers 
1 ^ k ^ 2. In the standard Milky Way MC model presented 
so far, the saturation amplitude of the ID velocity dispersion 
is 2.4 km s" 1 . 

We test the effect of velocities by performing experi- 
ments on two different sets of models for which the veloc- 
ities differ from the fiducial n300 model. The first set of 
models have different driven turbulent velocities. In these 
experiments, the other parameters of the simulations are 
identical to that from model n300 - only the saturated-state 
(3D rms) velocities differ by factors of 0.2 and 2. The sec- 
ond set of models are simply the original n300 models, but 
for which the velocities in each zone are manually modi- 
fied. As discussed in Paper I, the chemical evolution is not 
strongly affected by the velocity field. It is the metallic- 
ity, density, and background UV radiation field which are 
the primary factors determining CO form ation (see also 
iGlover et al.1l20ld ; Idover fc Mac Lowll201ll ). Given the in- 
sensitivity of molecule formation to the velocity field, we can 
directly modify the velocities in model n300. We thus con- 
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Figure 13. The linewidth-size relationship from model n300. 
Structures are identified directly from the 3D simulation. The 
lincwidths a are computed by taking the dispersion of the (ID, 
or Vz) velocities from the identified structures in the simulations 
(black, blue, and green) or from a Gaussian distribution with a 
dispersion of 2.4 km s — 1 (red). 

sider models with chosen velocity dispersions, where all the 
other parameters are equivalent to those in model n300. 

4-4-1 Different levels of turbulence 

Figure [13] shows the linewidth-size relations for clumps in 
the MHD models. "Clumps" are identified directly from 
the sim ulation (in the 3D dens ity cube), through dendro- 
grams (jRosolowskv et all |200S| ). The dendrogra algo- 
rithm identifies contiguous structures through iso-density 
contours. The linewidth of each clump is computed from 
the dispersion in the corresponding region from the 3D v z - 
cubel^J As evident in Figure [TUl the linewidth-size relation- 
ship can be reasonably expressed as a power law a oc R a . 
For the original n300 model (black circles) the best fit ex- 
ponent a=0.45. This is in goo d agreement to the observed 
linewidth size relationsh ip (e.g. [Larson 1981; Sol omon et al.l 
1 19871 ; I Hever et al.ll2009h . At small radii approaching the res- 
olution limit of the simulation, the velocities approach a 
minimum threshold corresponding to the microturbulent ve- 
locity of 0.5 km s" 1 . This chosen value of the microturbulent 
velocity is greater than the thermal velocity due to the ob- 
served linewidths at ~0.1 pc scales. At large scales, the ID 
linewidths ~2.4 km s _1 are the overall dispersion in the LoS 
velocities. 

The blue and green points in Figure [13] show the a - 



R relationship from simulations with different forcing am- 
plitudes, with 3D -Urms~ 1 and 10 km s _1 , producing LoS 
velocity dispersions of 0.71 and 5.8 km s _1 , respectivelvF 2 ! 
For these simulations, the best fit power laws produce a = 
0.08 and 0.61, respectively. The red points in Figure [LSl show 
the linewidth-size relationship from a model similar to n300, 
but with the velocities drawn from a random distribution 
with ct=2.4 km s" 1 (which produces a=0). 

Figure [T4l shows the relationship between X, Nu 2 > an d 
iVco from the original n300 model, along with the models 
with different velocity fields. In Figure 114b - f, the veloci- 
ties vary due to a different turbulent forcing amplitude, or 
simply due to replacing the velocities from the n300 MHD 
simulation with a Gaussian distribution. As discussed, the 
original velocities in model n300, (Fig. 114b .) have a LoS dis- 
persion of (7=2.4 km s _1 , due to turbulent forcing with 3D 
(rms) mean dispersion v Iuls — 5 km s _1 . 

Figure I14f c) and (e) show the X factor from models 
where the turbulent forcing is varied from the original n300 
to yield am ~ 0.71 km s , or am ~ 5.8 km s _1 , respec- 
tively. Figure [l3]demonstrates that the velocities of the large 
scale structures in these model vary by a factor of ~3 above 
and below the original n300 model, but are rather similar 
on the smallest scales due to the imposed microturbulence. 
Because the turbulent driving varies, the gas in these models 
contain different amounts of CO, and different temperatures. 
Comparing the other properties of these models (in Table[T|, 
the lower velocities in the 0.71 km s" 1 model result in lower 
temperatures ({T) vo i = 24 K), and lower amounts of total 
CO (as well as H 2 , (n H2 )voi = 127 cm" 3 ). The 5.8 km s" 1 
model, on the other hand, has (T) mass = 27 K. Neverthe- 
less, the X factor for these models are both within ~ 50% of 
the n300 value ((X) = 3.4 x 10 20 , 2.2 x 10 20 , and 1.6 x 10 20 
cm -2 K _1 km -1 s, respectively, for am = 0.71 km s -1 , 2.4 
km s~\ and 5.8 km s -1 .). 

4-4-2 Purely Gaussian Velocities 

The right panels in Figure [T3] shows the X factor from 
model n300, but with the velocities replaced with random 
values drawn from a Gaussian distribution with a— b) 2.4, 
d) 24, and e) 240 km s~\ For the a=2.4 km s" 1 model, 
the linewidth-size relationship differs significantly (Fig. I13|l 
from the fiducial model n300. The dispersions are similar 
on large scales, by design. However, at smaller scales the 
dispersions in the original model decreases, whereas there is 
a constant dispersion on all scales in the modified model. 
Nevertheless, the X factor distribution is very similar to the 
original model n300. The globally averaged X factor from 
the original model is recovered. That the X factor is very 
similar to the original model suggests that the details of the 
velocity structure, and its relationship to other physical prop- 
erties such as mass or size, do not play an important role 
in determining the X factor. In particular, clouds with very 
different linewidth-size relationships from a oc J? 1//2 may 
produce A~XGai- 

In models with random velocities having ID dispersions 



10 https:/ /people. ok. ubc.ca/erosolo/dendro. html 12 These LoS velocity dispersions include a contribution from the 

11 We use the z-component since the CO emission is calculated 0.5 km s — 1 microturbulent velocity, so are not exactly equal to 



for the i-direction. /\/3. 
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Figure 14. The variation of the X factor with Nn 2 and NqOi as m Fig. [8] but for models with different internal velocity fields, a) 
Original model n300, which has u\£>=2A km s —1 , b) model n300, but with velocities replaced from a Gaussian distribution with <r=2.4 
km s —1 , c) model with ctijj = 0.71 km s -1 (« r ms=l km s — 1 ), d) n300, but with Gaussian velocities with ct=24 km s — 1 , e) model with 
""in = 5.8 km s — 1 (« r ms=10 km s _1 ), and f) n300, but with Gaussian velocities with <r=240 km s _1 . In b)-e), the X factor distribution 
from the fiducial n300 model in a) is shown as the thin dashed histogram. The emission weighted X factor {X) is indicated in each panel. 



(7=24 and 240 km s 1 , the X factor is systematically low- 
ered, to (X) = 6.4 x 10 19 and 4 x 10 19 cm" 2 KT 1 km" 1 s, 
respectively. Therefore, the integrated CO intensity W is not 
sensitive to the velocity structure, but only the extent of the 
range in velocities. This occurs because with larger velocity 
differences between regions along a LoS, more CO line pho- 
tons are able to escape the cloud and ultimately be detected, 
resulting in an increase in W . Increasing Av thereby effec- 
tively reduces the percentage of mass that has r > 1 (and 
is therefore "invisible"), so that more CO becomes visible. 
Notice that the difference in the X factors of the models 
with dispersions of 24 and 240 km s _1 (FiglTlh and f) is 
modest. Thus, X does not simply scale inversely with a. In 
fact, the results from models with a — 2.4, 5.8, and 24 km 
s _1 show a behavior closer to X oc a" 1 / 2 than X oc cr _1 . We 
return to this point in Section f4. 61 

Of course, the range in velocities must be sufficiently 
well-sampled, so that there are no significant velocity gaps, 
which will indeed be the case in MCs. The cr=24 and 240 
km s _1 models have a large range in velocities, and con- 
sequently, likely have under-resolved velocity gradients. In 
the Appendix, we discuss how radiative transfer calculations 
may provide inaccurate intensities in regions where the ve- 



locity gradients are poorly resolved, and how our analysis 
accounts for this effect. 



4.5 Does "cloud counting" result in a constant X 
factor? 

One explanation for the lack of variation in the X factor in 
the Galaxy is that the integrated CO intensity is a measure 
of the number of optically thick "cloudlets" along the LoS, 
he reafter LoS. Th i s is k nown as the "mist" model proposed 
bv lSolomon et al l j 19871 ). Since we have knowledge of all the 
relevant quantities affecting the X factor, we can test this 
idea by inspecting the characteristics of the gas contributing 
to the observed emission. 

Figure [151 shows the spectrum from two LoSs. The top 
panels show Tb and r at each spectral channel. The bottom 
panels show the corresponding volume density and velocity. 
For reference, the observer is situated beyond LoS position 
in panels c-d, so that observed line photons are traveling 
from high LoS position towards lower LoS position (to the 
left on the plots). 

Both LoSs A and B have similar total column densities, 
1.0 and 1.1 xlO 20 cm -2 , respectively. However, the inte- 
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grated intensities differ by a factor of ~ 2.5. Perhaps sur- 
prisingly, LoS A, which has slightly lower column density, 
has a much higher intensity. As a result of this difference in 
W, the resulting X factors for these LoSs as computed from 
Equation [1] also differ. 

The differences in the total intensities and line shapes 
from the LoSs depicted in Figures I15fa -b can be under- 
stood by inspecting the density and LoS velocity distribution 
shown in Figures [Tob -d. LoS A has relatively low density gas 
£ 100 cm" 3 in the majority of positions along the LoS, at 
position 0-17 pc. Near LoS position 17 pc, there is a sharp 
jump in density, and an associated perturbation in velocity. 
The lowest LoS velocities associated with this shock, ,6—5 
km s _1 , are unique to this particular region along the LoS 
- no other region along this LoS has similar velocities. Ac- 
cordingly, observed emission in the —5 to —7 km s" 1 range 
originates from this high density, optically thick cloudlet. 
Emission at velocities —3 to 4 km s _1 from this cloudlet 
may be absorbed by the gas lying along the LoS with simi- 
lar velocities. However, this gas has very low column density, 
and so in fact does not significantly attenuate the emission 
from the shock. Thus, the whole line profile in the —7 to 
3.5 km s _1 range is due to the high density shock. Only the 
weak remaining emission at velocities Si 4 km s _1 is due to 
gas at positions ,$ 3 pc. We discuss some numerical effects 
in regions of such high velocity contrast in the Appendix. 

The LoS in Figures 115b does not have any gas with 
71h 2 > 2000 cm -3 . There are numerous optically thick re- 
gions, with a significant overlap in velocities. Some regions, 
such as the peak near position 8 pc, have small velocity gra- 
dients, while others contain a large range in velocities, as 
the peak near position 14 pc. Due to the large overlap in 
velocities, especially at —4 km s" 1 , the integrated optical 
depth reaches very high values up to 170. The combination 
of velocity overlap and lower density gas results in LoS B 
having a lower integrated intensity than LoS A. 

Figure [16] shows two other LoSs with lower column den- 
sities; these have equivalent Nu 2 = 3.5xl0 21 cm" 2 . Yet, the 
integrated intensity varies by a factor of ~3, resulting in 
an equivalent discrepancy in the X factor. Judging by the 
optical depth profile, there appears to be either one or two 
cloudlets. However, the detailed velocity and density profiles 
show that the structure is much more complex. 

Almost all of the gas along LoS C has densities Si 10 
cm -3 . Most of the gas velocities, especially those associated 
with the density peaks with iih 2 ~ 100 cm -3 , lie in the 
range —3 to —1 km s _1 . There are many density peaks with 
overlapping velocities, but the resulting line profile has only 
three peaks. Since r > 1, the observed intensity at a given 
velocity emerges from the last density peak along the LoS 
with the given velocity. Thus, much of the emission from 
LoS position > 10 pc is absorbed. By rerunning the radiative 
transfer on this LoS by excluding some of the gas, we have 
verified that a significant number of CO line photons are 
absorbed. We provide examples demonstrating this scenario 
in Section [4~6l 

For LoS D, a fraction of the gas has nn 2 ~ 10 cm -3 . 
This LoS has two well defined peaks in the line profile. These 
peaks are centered on —3 and 3 km s - . At those velocities, 
there are two distinct cloudlets with riH 2 ~ 1000 cm -3 . 
Clearly, the emission from LoS D comes from the most dense 
regions along the LoS. Since they are well separated in ve- 



locity, both are easily detected and contribute to a larger 
integrated intensity than LoS C. 

The comparison of observed profiles, densities, and ve- 
locities in Figures [151 and [T6l suggest that the simple "mist" 
model does not accurately capture the complexity intrinsic 
to line radiative transfer from a turbulent medium. In the 
comparison of LoS A and B in Figure [T5l both of which have 
a large Nu 2 , A has one true cloudlet with very high density. 
This cloudlet is only present in a localized region along the 
LoS, but contains a large range in velocities, and therefore is 
the source of most of the observed emission. Other emitting 
regions with different velocities only contribute slightly to 
the spectrum, due to their significantly lower densities. LoS 
B has numerous lower density cloudlets, many of which have 
similar velocities. The integrated intensity is lower, though 
the total amount of gas along this LoS is slightly larger than 
LoS A. 

In the comparison of low density LoSs in Figure ITfJl both 
LoSs have a few regions with clear density peaks. However, 
the dominant cloudlets in LoS D have larger density, and 
span a larger range of velocities. Thus, LoS D has a higher 
total intensity than LoS C even though the total column 
density in both LoSs are equivalent. 

The "mist" model would predict that LoSs with more 
high density cloudlets would have higher integrated intensi- 
ties, since the cloudlets are separated in velocity. However, 
comparison of LoS A and B show the opposite: one very 
dense cloudlet with a large velocity gradient may be respon- 
sible for the whole LoS profile. This integrated intensity W 
may be larger than that from a different LoS with a larger 
number of cloudlets, but nevertheless a similar total column 
density and total velocity width. Further, two LoSs with 
equivalent Nh 2 an d similar number of cloudlets may have 
different W if one LoS has overlapping cloudlets in velocity 
space, while the other LoS has well distributed cloudlets in 
velocity space. This effect is partly at work in the compari- 
son of LoS C and D in Figure 1161 

We have shown that in a turbulent medium, emit- 
ting cloudlets do in fact overlap in ve l ocity space (see als o 
iBallesteros-Paredes fc Mac Low! |2002| ; IShettv et ail l201Ch . 
and that individual cloudlets may have a range in veloc- 
ities, thereby dominating the emission at all velocities in 
that LoS. These factors all compromise any simple relation 
between W and Nn 2 , and therefore any direct scaling be- 
tween X and Nh 2 ■ Of course, this analysis considers indi- 
vidual LoSs through a molecular cloud, whereas the idea of 
a constant X factor is usually discussed in the context of 
whole molecular clouds. Accordingly, we next consider the 
cloud-averaged CO spectrum. 

4.6 Cloud-averaged X factor and spectra 

Figure [171 shows the averaged spectrum for the whole n300 
cloud model, as well as model n300 where the velocities are 
are manually replaced with a = 24 km s _1 (see Fig. 1141 0. 
Each point shows the mean brightness temperature of each 
channel in the synthetic observation. The integrated inten- 
sity of this spectrum is W— 40.2 and 137.5 K km s _1 for 
the fiducial and velocity-altered model, respectively. 

The best fit Gaussians to the averaged spectra are also 
shown in Figure [171 The spectra are very well fit by Gaus- 
sians. Using the peak and dispersion of the Gaussian fit, the 
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Figure 15. Top panels - CO spectra (lines - left axis)) and optical depth (circles - right axis), and bottom panels - H2 column densities 
(thin - left axis) and LoS velocities (thick - right axis) from two LoSs in the n300 simulation. The corresponding integrated intensities 
and X factors are listed in a-b, and the column densities are listed in c-d. 
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Figure 16. Similar to Fig. 1151 but for different LoSs. The column density of both LoSs are equivalent. 



corresponding integrated intensities reproduce the true in- 
tegrated intensities. Using W along with the mean column 
density N h 2 = 9.0 x 10 21 cm -2 in Equation [T] results in 
X= 2.2 xlO 20 and 6.5xl0 19 cm" 2 K" 1 km" 1 s for the orig- 
inal and 24 km s _1 models, respectively. These values of the 
X factor are equivalent to the emission-weighted mean X 
found in Section 4.1 and listed in Figure [T4l The X factor 
value for n300 is nearly identical to the reference value in Ta- 



ble [T] For the altered velocity model, {vl) ma ss = 24 km s 1 

so the reference (X) re f = {n H2 }voiL/[{T} m ass({v^) m asa)^] 
= 1.9 xlO 19 cm -2 K" 1 km" 1 s. This value is significantly 
lower than X determined by the spectrum, largely due to the 
fact that Tb / {T) m ass is a factor 3 lower for the 24 km s _1 
model than for model n300 with <jid=2.4 km s . Therefore 
the mass-weighted quantities {T} maS s and (f 2 )m 2 S3 , cannot 
fully account for W . We note that the line widths of the 
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Figure 17. Averaged spectrum from model n300 (black), and 
model n300 for which the velocities were replaced from a random 
distribution with a = 24 km s _1 (red). Lines show best fit Gaus- 
sians. Integrated intensity, fit Gaussian parameters, and X factor, 
are listed for both models. 



Figure 18. Averaged spectrum from original model n300 (black 
circles), along with spectrum from the original model but only 
gas with nn 2 < 1000 cm~ 3 (red squares) and > 1000 cm" 3 (bluc 
triangles) is considered. Lines show best fit Gaussians. 



observed spectrum are slightly larger than the intrinsic ve- 
locity dispersion. This is due to some extent to our choice 
of the microturbulence, which we discuss in more detail in 
the Appendix. 

The cloud- averaged spectra in Figure [T7] show very 
smooth profiles, and do not demonstrate any of the vari- 
ability seen along individual LoSs. Clearly, the cloud with 
a larger velocity dispersion is brighter because the emission 
from more gas with a larger dispersion is able to escape the 
cloud. Even though both models have the same mass, their 
integrated emission differs by a factor 3.4. This is due to 
a combination of a larger velocity dispersion (by a factor 
~ 10), and a lower peak brightness temperature (by a factor 
~ 3) for the 24 km s _1 model compared to model n300. As a 
result, use of the integrated intensity alone, or any X factor, 
would provide inaccurate mass estimates. 

We have demonstrated that there is no simple relation- 
ship between W and Nu 2 in highly turbulent clouds (see also 
Paper I). This lack of correlation is due to the combination of 
the optically thick nature of CO, and the complexity in the 
structure of CO bright regions: emitting gas along individ- 
ual LoSs sometimes overlaps in velocity space and sometimes 
does not, and even lower density gas not necessarily found 
in cloudlets contributes significantly to the observed inten- 
sity. Under such circumstances, there should be significant 
amounts of gas that is untraceable in the CO line. 

Figure [18] shows the averaged spectrum of three models. 
The black spectrum is from the original n300 model. The 
red and blue points show the results of the radiative transfer 
calculations from model n300, but only when gas with 71h 2 < 
1000 cm -3 and > 1000 cm -3 is included in the calculation, 



respectively. The CO and H2 density is set to zero in zones 
with nH 2 < 1000 cm~ 3 or > 1000 cm" 3 , respectively. 

The gas mass in the model that has hh 2 < 1000 cm" 3 is 
5.3 x 10 4 Mq , which is roughly 2/3rds of the total mass in the 
original n300 model. However, the integrated intensity from 
this model is 90% of that from model n300. The total mass at 
7iH 2 > 1000 cm" 3 , 2.5xl0 4 Mq, accounts for the remaining 
1/3 of the mass in model n300. But the integrated intensity 
from this model is 26.2 K km s , which is 65% of the total 
W of the original model. 

Due to the significant amount of self-absorption occur- 
ring in both low density and high density gas, removing 
some gas simply allows other gas along the LoS to become 
"visible." There are LoSs where some CO emitting regions 
go completely undetected, since other optically thick gas ly- 
ing in the same LoS with similar velocities will absorb its 
emission. Consequently, a spectrum from a given LoS may 
depend on the observing direction. For example, if the ob- 
server were located at LoS position >20 pc in Figures [15] and 
1161 instead of at LoS position <0 pc, the spectra may have 
different profile shapes. We have indeed found such discrep- 
ancies in many LoSs. This is further confirmation that there 
is no simple correlation between A^and W . We stress again 
that the use of a constant X factor for numerous clouds will 
not provide reliable estimates of the molecular mass. 
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5 DISCUSSION 

5.1 Why is the Galactic X factor found to be 
nearly constant? 

After investigating the dependence of the X factor on the 
physical characteristics of synthetic MCs, we are now in a 
position to interpret the observed trends, as well as compare 
our results with previous efforts. 

For our model with Milky Way parameters, we have 
demonstrated that under certain conditions, using the mass- 
weighted mean values of temperature and velocity dispersion 
may provide reasonably good estimates of X factor. Thus, 
it might seem unsurprising that systems with similar mass- 
weighted temperatures or velocity dispersions are able to re- 
produce the X factor from the original model. By exploring 
models with different temperatures and velocities, however, 
we find that X re f = (Nn 2 ) /({T) mass a v j os ) does not in gen- 
eral agree with the value of (X) =Nu 2 /W computed from 
radiative transfer models. Further, (X) does not simply de- 
pend on the inverse of the temperature or velocity. 

Even though the experiments testing the variation of 
a constant kinetic temperature on the CO line do not self- 
consistently follow the coupling of other MC properties with 
T, we may compare our results with previous investigations. 
We find that X oc T -0,5 , which is a weaker X—T dependence 
than the X oc J 1-1 - 3 behavior found by iKutner &: Leund 
(1985). Their models have velocities which are dominated 
by microturbulence, but nevertheless obey the empirical 
linewidth-size relationship. We find that microturbulence 
does not significantly affect the integrated intensity, nor cor- 
respondingly the X factor (see Appendix showing profiles 
with different microturbulent velocities). In our models, the 
CO intensity is more sensitive to the macroturbulent veloc- 
ities. 

Thus, our results are in cl oser agreement to the investi- 
gation by IWolfire et all (|l993l ), who find that macroturbu- 
lent or clumpy models are better able to reproduce observed 
line profiles than microturbulent model s. However, a key dif- 
ference between our model an d that of lWolfire et all (l993, 
as well as lDickman et akll 19861 ) is that in their clumpy mod- 
els, much of the CO emission emerges from distinct optically 
thick clumps. Further, their models require approximately 
one clump along each LoS at a given velocity in order to 
reproduc e observed sp e ctra. 

The IWolfire et all (|l993h models are conceptually sim- 
ilar to the "mist " model often employed to explain the con - 
stant X factor (|Dickman et al.l 1 19861 ; ISolomon et al.1 Il987l ) 
and described in Section 4.5. In our turbulent models, only 
the cloud-averaged CO linewidths provide accurate mea- 
sures of the total (ID) velocity dispersion. Along an indi- 
vidual LoS, the line shape is strongly dependent on the den- 
sity and velocity structure. We find that though dense gas is 
responsible for much of the observed CO emission globally, 
along individual LoSs diffuse gas may also be a significant 
source. Depending on the gas velocities along the LoS, lower 
density regions may also be optically thick, and so may con- 
tribute substantially to the emergent intensity. 

One of the key results of our work is that the X factor 
is insensitive to the detailed linewidth — size relationship in 
the cloud. In Section 4.4 we demonstrated that a cloud con- 
sisting of gas with random velocities has an almost identical 
X factor distribution with a cloud that has a well defined 



a oc linewidth-size relationship. This suggests that the 
linewidth-size relationship has very little influence in deter- 
mining the X factor, or W. The lack of dependence of W 
on the detailed velocity structure suggests the use of a rep- 
resentative velocity width for the J dv term in Equation 
[4] Although linewidths computed from the mean spectrum 
may differ slightly from (vl)\lass, a more significant effect of 
radiative transfer is to reduce Tb when the velocity width 
increases. Consequently, we find that W does not increase 
linearly with the velocity dispersion. In particular, we find 
that an increase of the velocity width by a factor 10 or 100 
only reduces (X) by a factor 3 or 5 relative to model n300 
(see caveats in Appendix). 

Along with Nh 2 , W determines the X factor through 
its traditional definition given by Equation [1] Based on our 
tests, for A r H 2 ~ 10 22 cm -2 , the X factor can only vary by 
a factor 5, given the limited range of dense-gas temper- 
atures 10 - 25 K a nd velocities ^2 - 5 km s" 1 observed for 
Gala ctic MCs (e.g. lSolomon et al.lll987l ; lRoman-Duval et al.l 
l2010h . The mass-size relationship proposed bv lLarsonl (|l98ll ) 
implies a constant mean column density in MCs. Recent 
observations have challenged the notion that all MCs have 
column densities, averaged over large scales, similar to one 
another (|Hever et al.ll2009l ; iKauffmann et~al] l2010bl ). Nev- 
ertheless, column densities are not expected to vary signif- 
icantly from cloud to cloud, and they generally fall in the 
range from 10 21 -10 22 cm . Structures with much lower col- 
umn densities than ~ 10 21 cm' 2 probably would not be ob- 
served as molecular clouds, as there would be insufficient CO 
to be detected. Clouds with column densities much larger 
that ~ 10 22 cm -2 would be very gravitationally unstable 
- rapid collapse, fragmentation, and the destructive effects 
of star formation would all prevent clouds from maintaining 
such high column densities over long periods of time (within 
an environment of much lower mean density and column 
density). Thus, we believe the limited range in Nu 2 , velocity 
dispersion, and MC temperature together explain the lim- 
ited range in the X factor found by numerous observations 
of Galactic MCs. 



5.2 Does a constant X factor require virialized 
clouds? 

A common interpretation of the limited range in X factors 
from observations of Galactic MCs i s that they are "viri- 
alized" objects. The "mist" model of ISolomon et al.1 (|l987l ) 
invokes the "virialized" cloud depiction to explain the lim- 
ited range in the measured virial mass - CO luminosity 
(Mvt I Leo) ratio. Generally, a mass-to-luminosity ratio, a, 
is simply the X factor written in terms of a gas mass and CO 
luminosity, instead of Nh 2 and W, so that X= 4.5 x 10 19 a 
if a has the units M©/(K km s _1 pc -2 ). Since the deriva- 
tion of a constant olvt =Mvt / Leo for v irialized clouds 
has been presented in numerous works (e.g ISolomon et al.l 
ll987l ; lMalonevlll990l ; lYoung fc Scovillelll99ll ). we only briefly 
outline these arguments, and then discuss how our models 
relate to such a description. 

For a cloud in virial equilibrium, the internal veloci- 
ties are fully governed by the gravitational field due to the 
gaseous mass, 

GM VT 



2 

°~VT 



511 



(8) 
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Here, we have used "VT" to label the one-dimensional virial 
velocity width ovt , and the factor 5 assu mes a spheri- 
cal, uniform cloud (jBertoldi fc McKedfl99^ . Note that in 
this notation, ovt is the ID velocity dispersion (assuming 
an isotropic distribution), whereas a refers to the observed 
linewidth along the LoS. 

Cloud-averaged spectra often have Gaussian shapes, so 
that the CO luminosity can be computed using the Gaussian 
line parameters Ts,o, the peak brightness temperature, and 
linewidth a, along with the projected area irR 2 : 



Leo = V2nT B ,ovnR ■ (9) 

Taking Equations [H] and together, and assuming a — avr, 

_ Mvr o/R 

a VT = -j oc - — . (10) 

Leo Jb,o 

Using Equation!] a/R oc (GM/R 3 ) 1/2 so that Equation [TOl 
is equivalent to 

1/2 

(11) 



_ M V T P 

avr = -j oc — — 

Leo Jb,o 



If p and Tb.o take on a limited range of values, this yields 
an approximately constant avr and X. 

When further combined with an empirical linewidth- 
size relation a oc R 1 ^ 2 , equation [8] yields 



Mvt oc a 



(12) 



and equation [9] yields Leo oc Tb,qcf 5 . These relations com- 
bine to give 



Mvt oc 



(13) 



CO observations are often used to obtain cloud masses 
through Equation [5] assuming a = ovt- Such an analy- 
sis usually produces a strong correlation between Mvt and 
Leo, with a power law index 0.8. Though we have only fol- 
lowed proportionalities in Equations [TU] - 1131 it is straight- 
forward to consider the coefficients e xplicitly (e.g. iMalonevi 
1199(1 ISolomon fc Vanden BoutfeOoEl) . The observed M VT - 
Lco trend also results in good agreement in these coeffi- 
cients. The agreements in both the power law exponents 
and coefficients are often taken as evidence that MCs are 
virialized objects, and that Nn 2 is approximately constant 
for MCs, since Nr 2 oc E oc M/R 2 ~ constant if a 2 oc M/R 
and a oc R 1/2 . 

The apparent agreement between the observed 
Mvt I Leo and the analytical derivation outlined above de- 
pendends on two crucial assumptions, the mass-size and 
linewidth-size relations. As discussed above, numerous ob- 
servations demonstrate a a oc R 1 ^ 2 re l ationship over 
a ran ge of scales fo r MCs l|Larsonl Il98ll ; ISolomon et al.l 
19871; iBolatto et al.l l2008h . However, recent analysis by 
Hever et al.l 1)20091 ) of higher resolu t ion o bservations of the 
clouds studied by ISolomon et al.l l|l987h found that the 
coefficient of the linewidth-size relationship may not be 
constant. Furthermore, the observed mass-size relationship 
M oc R 2 is rather uncertain. Using observations of MCs 
in 12 CO and 13 CO, which should trac e more of the cloud 
;as d ue to the lower optical depth, iRoman-Duval et al.l 
2010l ) find an exponent of 2.4 rather than 2. By assess- 
ing the mass on a con tinuous range of s cales in a num- 
ber of Galactic MCs, iKauffmann et all (|2010al lbh found 



no universal power law mass-size relationship. Addition- 
ally, when optically thin (e.g. dust) observations are em- 
ployed, projection effects may lead to overestimated masses 
teallesteros-Paredes fc Mac Lowll2002l ; iGammie et~al]|2003l ; 
IShettv et alj|2010h . Another potential issue is that the defi- 
nition of cloud size R may affect any M estimates. Usually, 
the size R is set to the radius of a circle with area equivalent 
to that of the projected area of the cloud. Using a different 
definition of R may result in an altered mass-size relation, 
and possibly even a different linewidth-size relation. These 
caveats will all affect any derived mass - luminosity correla- 
tions. 

In this paper, we have studied CO emission from turbu- 
lent MC models which do not have self-gravity, and thus are 
not "virialized." Nevertheless, the virial mass inferred from 
the synthetic CO observation accurately recovers the cloud 
mass (7.8xl0 4 M Q ) to within 1%, for the n300 model. Since 
the parameters adopted for the n300 model were designed to 
represent Milky Way MCs, this agreement may lead to the 
spurious inference that the MC is "virialized." Additionally, 
we have found that the integrated emission, and hence the 
X factor for the adopted value of Nn 2 , is insensitive to the 
velocity structure (see Section 4.4). As a result, a cloud with 
a S> ovt or <7 < o vt which is far from "virialized" could 
still have quite similar X values to a cloud with a ~ ovt- 
Co nsequently, our analy sis suppo r ts the arguments put forth 
by IMalonevi l| 19901 ') and ICombesI ([1991) that a constant X 
factor does not require MCs to be "virialized." 

In effect, our analysis shows that X^XgbI is simply a 
result of the limited range of both the numerator and denom- 
inator of Equation [4] or of the parameters entering the right 
hand side of Equation [10] or [11] for gravitationally-bound 
MCs: the densities, brightness temperatures, and (observed) 
velocities. Only if the velocities are governed solely by gas 
self-gravity through Equation and a universal linewidth- 
size a oc R 1 / 2 relationship holds can a (nearly) constant 
Mvt I Leo ratio be expected for all MCs. The linewidths of 
observed clouds, which measure the internal motions, may 
or may not be set solely by gas self-gravity. An additional 
caveat is that the Equation [8] may not be sufficient for in- 
demnifying the dynamical state of a cloud, since it does not 
account for surf ace, magnetic and time-va rying terms in the 
virial equation ( Ba llesteros-Pa rcdcs 20061) . 

5.3 Future work: X factor in molecular dominated 
regions 

For molecular environments different from Milky Way 
GMCs, the different (volume and column) densities, veloci- 
ties, and/or temperatures, should contribute to a variation 
in the corresponding X factors. As indicated in the previous 
section, we have focused on one model with densities similar 
to those of Milky Way MCs. Molecular dominated regions 
such as ULIRGs and the Galactic center have much higher 
column densities, in addition to elevated temperatures and 
observed velocity dispersions, so the results found in this 
work are not directly applicable to those sources. In this 
section, we briefly explore the results from the previous sec- 
tion in the context of ULIRGs, and motivate an analogous 
study of higher column density models. 

In molecular dominated regions, t he X factor is believe d 
to be lower than XGal (see e.g. Fig. [T1 iTacconi et al.ll2008h . 
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In part, this is because the standard value of Xcai yields 
molecular masses that are sometimes larger than the dy- 
namical masses of the galaxies. In orde r to avoid this un- 
physical situation, iDownes et al.l (|l993l ) proposed that the 
masses traced by CO observations of ULIRGs are not the 
virial masses, but rather the geometric mean of the gas and 
dynamical masses, i.e. Leo oc T B ,o{M g asM d y n / p ga s) 1/2 ■ In 
their interpretation, the CO linewidths trace the dynamical 
masses of the starbursting, molecular rich systems. Thus, 
the virial velocities and masses in Equation [8] must be re- 
placed by the dynamical velocities and masses of ULIRGs. 
Instead of using avrLco (Eqn. Hip to estimate masses, 
one can estimate the product of dynamical and gas mass of 
ULI RGs through (aLgo) 2 , where a = 2.6n],L 2 /TB,o (Eqn. 
5 in lSolomon fc Vanden Boutll2005l ). The ratio M gas /L C o 
computed from ULIRGs is ~ 3 — 6 times lower than th at 
estimated from MCs (e.g. ISolomon fc Vanden Bout| [2005). 

ULIRGs have observed CO li newidths of the or- 
der of a few hundred k m s _1 (e.g. ISolomon et al.l 1 19971 ; 
IDownes fc Solomon! |l99Sh . Since the molecular component 
of ULIRGs are inclined rotating disks, the observed CO 
linewidths are primarily tracin g the rotational mo t ion o f 
the molecular-dominated ISM. IDownes fc Solomon! (1 19981 ) 
find that turbulent velocities of ULIRGs must be signifi- 
cantly higher than Galactic MCs in order to accurately fit 
their observed lines. However, due to the compact extent 
of ULIRGs and the high optical depth of CO, accurately 
separating the turbulent and rotational contributions to the 
observed linewidths is not trivial. 

The observed brightness temperatures of ULIRGs are 
significantly lower than Galactic MC values. However, 
ULIRG kinetic temperatures can be £ 1 times larger 
due t o enhanced star formation activity (e.g. lSolomon et al.l 
1 19971 ). In Section 14.61 we showed that the brightness tem- 
perature indeed decreases with increasing velocity disper- 
sion. Accurately quantifying this Tb— a relationship calls for 
more realistic models which self-consistently contain higher 
temperature, higher dispersion, and have sufficiently high 
resolution (see Appendix). 

In Section 14.41 we showed that a model with random 
velocities with dispersion 240 km s~ , similar to the val- 
ues found from ULIRGS, produces X= 4.0xl0 19 cm" 2 K" 1 
km" 1 s. In addition, we found X= 7.9xl0 19 cm" 2 KT 1 
km -1 s for the T=100 K model from section 4.2. These 
X factor values are a factor ~3-5 lower than -Xcai and the 
value found from the original n300 model, and are consistent 
with X factor estimates from ULIRG observations. These 
could be due to either extremely high small-scale turbu- 
lent velocity dispersions (as experienced by photons escaping 
from the molecular disk) and moderate gas temperatures, or 
much higher gas temperatures and moderate small-scale tur- 
bulence combined with very large non-circular motions on 
10 2 — 10 3 pc scales (still within the typical beam size for 
molecular observations, but larger than the molecular disk 
thickness) . We note that none of our current models has col- 
umn densities comparable to ULIRGS. The effects of higher 
star formation, very high density and total column density, 
and rotation all need to be studied in detail to understand 
the CO-H2 conversion factor from ULIRG systems. More 
detailed modeling should quantify the contribution of ro- 
tational and turbulent velocities, along with other physical 
properties in setting the X factor of ULIRGs. Efforts con- 



sidering such environments are currently underway (current 
authors, and Narayanan et al. in prep). 

6 SUMMARY 

6.1 Overview: Interpreting the X factor 

By performing radiative transfer calculations on model MCs, 
we assess the physical properties which are most responsible 
for setting the X factor, the ratio of molecular column den- 
sity to integrated CO intensity. The fiducial model MC is 
a result of a magnetohydrodynamic simulation of a gaseous 
medium responding to driven turbulence with a treatment 
of chemistry which tracks the formation and destruction of 
H2 and CO; the resulting MC has many properties similar to 
those found in Galactic MCs (Section and Paper I). After 
discussing the X factor from the fiducial model, we modify 
the physical characteristics of the model MC and recompute 
the emerging CO line emission from the altered models. We 
then compare the modified X factor with that provided by 
the original model (Section|4}, and discuss our results in the 
context of observed trends (Section [5]). 

Our analysis was aimed at understanding the limited 
range in the X found from Galactic MC observations, shown 
in Figure [T] We discussed the trend of a higher X factor for 
low density regions in Paper I. Basically, environments with 
lower CO densities have lower CO intensities. This could be 
due to lower total density, higher background UV radiation 
fields, or lower metallicities. Though S gas in the simulations 
ranges between ~50-500 Mq pc -2 , and the X factor along 
individual LoSs can vary significantly, in all except the very 
low metallicity models the cloud- averaged (X) ~ XgoI (see 
Fig. [2J). In this work, we have found that it is the limited 
range in velocities and temperatures that constrains W, and 
thereby X, for a given column density. 

This interpretation applies for systems above a thresh- 
old column density beyond which CO is optically thick (~20 
Mq pc -2 ), but below very high column densities commonly 
found in ULIRGs. The precise value of the threshold col- 
umn density depends on a number of environmental factors 
of the source, such as the metallicity, UV radiation inten- 
sity, and dust density. Below this column density, CO is 
optically thin, so a small decrease in column density results 
in a stronger decrease in CO intensity W , and thereby an 
increase in X. At larger column densities, CO is optically 
thick. Thus, an increase in CO abundance may or may not 
lead to an increase in CO intensity, depending on the dis- 
tribution of gas velocity. Much larger CO column densities 
are usually found in molecular rich sources such as ULIRGs, 
which generally have higher temperatures, total densities, as 
well as observed velocity dispersions (due to rotation and/or 
turbulence) — conditions which are not self-consistently fol- 
lowed in the models presented here. 

An important consequence of our conclusion is that a 
given X factor value may not always provide an accurate 
mass estimate of MCs detected from CO observations. Ob- 
servationally, X factor relations between total mass and 
CO brightness may coincidentally work well over suitably 
large areas of molecular clouds. This may simply reflect 
the limited range in the mean column density, temperature, 
and velocity dispersion within typical MCs. These underly- 
ing global properties of the emitting regions make a bridge 
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between the average of many LoSs with individual LoSs, 
where the relationships generally fail. MC masses estimated 
through a chosen constant X factor should only be thought 
of as a first approximation of the gas mass, to within ~50%. 
Indeed, that derived X factors for Galactic MCs themselves 
vary by up to a factor of five has already been a warning 
sign that a universal CO luminosity to gas mass conver- 
sion does not exist (e.g. lSolomon et al1ll987l ; IOka et al.lll998l ; 
I Pineda et al.ll2008l . and Fig.[T]). Due to its high optical depth, 
there may be a significant amount of CO self-absorption. As 
a result, CO intensities only provide rough estimates of the 
gas mass. 



6.2 Main Results 

1) We find a roughly constant X throughout the model MC 
only when the CO intensity is integrated over all velocities. 
When intensities at each PPV location are considered along 
with the column density Nn 2 ,u associated with that PPV 
position, there is a clear trend of an increasing X factor 
with increasing column density. This occurs because in high 
density regions the CO line is saturated, so increasing the 
column density only leads to more self-absorption, resulting 
in an increase in the column density - CO intensity ratio 
(Section 4.1). 

2) The variation in X between different LoSs through an 
individual cloud is not due to the variation in temperature 
within the MC. When setting the temperature to the mass- 
weighted average value, the original X factor distribution is 
recovered. The cloud-averaged X factor is thus determined 
primarily by the densest regions in the MC. We find a weak 
dependence of X on T: roughly (X) oc T" ' 5 from 20 K to 
100 K if density and velocity structure are fixed (Section 
4.2). 

3) Even though the X factor is defined with a column 
density, it also depends on the volume density tih 2 . This is 
due in part to the saturation of the CO line at high densities, 
as indicated in 1) above. Increasing the mean density from 
no — 300 to 1000 cm _3 only decreases (X) by 15%, however 
(Section 4.3). 

4) Similarly, at high CO abundances X increases with 
-/Vh 2 due to line saturation. At very low abundances (/co ~ 
10 -6 ), the line is optically thin, so W is directly proportional 
to -/Vh 2 producing a constant X factor. Since ther e is a range 
in CO abundances in MCs (|Glover et all l20ld . Paper I), 
the X factor shows signatures of both high (saturation) and 
low opacities (Section 4.3). On average, the X factor is not 
strongly dependent on the /co if /co is in the observed 
range G 10~ 5 - 10~ 4 (Section 4.3). 

5) The fiducial model MC has a linewidth-size relation- 
ship a oc ii 1 / 2 , and produces X~Xc a i, similar to observa- 
tions. Nevertheless, the X factor from models without any 
linewidth-size relationship, but with a similar range in veloc- 
ities, reproduces the X factor trends in the fiducial model. 
Further, in models with a larger range in velocities, there 
is less radiation trapping, leading to larger CO intensities 
and correspondingly lower X factors. Thus, it is the range 
in velocities that determine the X factor. The details of the 
velocity structure do not influence X (Section 4.4). Since 
the CO brightness temperature also depends on a, the X 
factor does not decrease linearly with increasing a, but in- 



stead shows a variation closer to X oc <r~ 1,/2 for a = 2 — 20 
km s _1 (Section 4.6). 

6) A CO-bright LoS with a large linewidth may be due 
to a single high density region. Alternatively, a similar line 
profile may result from a LoS with numerous lower den- 
sity gaseous peaks spanning a range of velocities. Therefore, 
two LoSs may have similar integrated intensities, but vastly 
different morphologies. Due to this diversity in morphology 
in turbulent clouds and the optically thick nature of CO, 
CO emission does not trace H2 in a one-to-one fashion on 
individual LoSs. Th e "mist" model (|Dickman et al.l 1 19861 : 
ISolomon et aflll987| ) does not accurately capture the com- 
plexity of turbulence and radiative transfer of the CO line 
in MCs (Section 4.5). 

7) Similar to individual LoSs, cloud-averaged spectra 
do not represent the total molecular content in a one-to-one 
fashion. Due to the blending of numerous emitting features, 
the averaged CO line is Gaussian. However, in experiments 
where some gas is removed, the intensity of the resulting 
spectra is not simply scaled down by the factor that would 
be appropriate if CO were a linear tracer of molecular gas 
(Section 4.6). 

8) Since we find that the X factor is insensitive to the 
detailed velocity structure, we argue that the reason Galac- 
tic MC observations generally find X~ Ao a i is due to the 
limited range in observed linewidths ~ 2-5 km s _1 , as well as 
in the temperatures within dense, CO-abundant regions in 
MCs (~10-25 K, Section 5.1), and in the column densities. 
Thus, the nearly constant X factor found in Galactic MCs 
does not necessarily require that a oc R 1 ^ 2 , or that clouds 
are "virialized." Although a near-constant Mvt / Leo ratio 
would result if the velocities were governed solely by gas self- 
gravity, and a oc R 1 ^ 2 universally, these are not necessary 
conditions for X to be nearly constant, and may or may not 
be true for Galactic MCs (Section 5.2). 

9) Given our explanation that the approximately con- 
stant X~Xoai found in Galactic MCs is due to the lim- 
ited range in MC properties, masses deduced through a 
constant X should only be considered as a rough esti- 
mate, within a factor ~ 2. Taking our results from Paper 
I and the current work together, we find cloud- averaged 
X ~ 2 - 4x 10 20 cnT 2 KT 1 km -1 s for Solar-metallicity mod- 
els with T ~ 10 — 20 K, velocity dispersions a ~ 1 — 6 km 
s _1 , and = Nh 2 ~ 2 — 20x 10 21 cm -2 - i.e. similar to Galactic 
MCs. The larger X factor values found in lower metallicity 
systems can be understood due to lower CO abundances 
(see references in Section 5.4). In molecular-dominated re- 
gions such as ULIRGs, the combined effect of higher CO and 
H2 volume and column densities, higher temperatures (due 
to enhanced star formation), and larger velocities (due to 
rotation and/or turbulence) all need to be self-consistently 
treated to assess the X factor. 
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Figure Al. CO spectra from LoS B (see Fig. ll5|l . Three different 
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and 0.75 (blue) km s — 1 . The velocity-integrated intensities W 
are indicated. 
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APPENDIX A: THE EFFECT OF 
MICROTURBULENCE 

As indicated in Section [3] for the radiative transfer calcu- 
lation, we set the microturbulent velocity t> m t r b= 0.5 km 
s" 1 in each zone of the simulation. This value is motivated 
by loosel y extrapolating the observed linew idth-size rela- 
tionship jLarso nf ll98ll ; ISolomon et ail 1 19871 ) down to the 
physical resolution of the simulation (~0.1 pc). Due to ob- 
servational uncertainties and the scatter in the observed 
linewidth-size relationship, the turbulent velocity at any 
given scale is not precisely known. We have thus explored 
how other values of u m trb may affect our results. 

Figure IA1I shows the spectrum from one LoS. It cor- 
responds to LoS B shown in Figure 1151 Three spectra are 
shown, with different values of v mt rb~ 0.25, 0.5, and 0.75 km 



s . With higher microturbulence the line is broadened, and 
gaps or sharp features in the spectrum can become smooth, 
as evident near velocities of —0.5 and —5.0 km s" 1 . 

The variations in the spectra lead to slight differences 
in velocity integrated intensity W. For w mtrb = 0.25, 0.5, and 
0.75 km s" 1 , W = 58.9, 64.2, and 73.0 K km s _1 . The dif- 
ferences in integrated intensity will of course result in cor- 
responding differences in the X factor. For this LoS, the 
column density is l.lxlO 22 cm -2 , providing X — 1.9, 1.7, 
and 1.5 xlO 20 cm -2 K _1 km -1 s. In the cloud- averaged 
spectra (e.g. Fig. I17p . the smooth Gaussian profiles are not 
significantly affected, and the linewidths only vary by ^ 
6% from the fiducial model (e.g. Fig. I17p . For w mtrb = 0.25, 
0.5, and 0.75, the peak Tb,o= 5.1, 5.8, and 6.3 K, leading 
to cloud-averaged X factors 2.6, 2.2, 1.9xl0 20 cm" 2 K" 1 
km" 1 s, respectively. 

We have found that the precise value of v m trb does not 
affect the cloud-averaged spectrum — and the X factor — 
significantly. The most noticeable difference occurs in the 
line shapes along individual LoSs. Nevertheless, even then 
the differences are minimal. Any value of f mtr b within the 
reasonable range G 0.25 — 0.75 would produce very similar 
CO intensities. Thus, the values of the X factor found in this 
work, and the conclusions drawn from its analysis, would not 
be affected if different values of v mtI b were employed. 



APPENDIX B: NUMERICAL RESOLUTION 

While performing radiative transfer calculations, there is a 
potential for obtaining inaccurate line intensities due to un- 
resolved velocity structure. In this section, we discuss why 
resolving the velocity gradients is important, and how we 
have accounted for under-resolved velocities in our study. 

In a grid-based model, the physical properties are only 
denned in each zone - either at the zone center or zone edge. 
While integrating the equation of radiative transfer, or ray- 
tracing, through a given LoS, the optical depth is only up- 
dated when stepping from one grid zone to the next. For 
a given "observed" frequency v, corresponding to velocity 
w bs, the optical depth r„ is only modified in zones with 
LoS velocities that correspond to the non-zero portion of 
the Gaussian profile function cj> v (see Eqn. [7] and Eqn. 6 in 
Paper I). 

As an example, consider the velocity profile of a sample 
LoS shown in Figure iBll The LoS component of the veloc- 
ity v s is only defined at discrete locations along the LoS s 
(e.g. edges of grid zones). At the indicated v h s , due to the 
large velocity gradients, there would be no emission (or, if 
optically thick, absorption) from this region, since none of 
the velocities overlap with Dobs- Of course, if the (Gaussian) 
width of 4>v were large enough (e.g. due to microturbulence), 
then there may be a contribution from one or more zones. 
For our present discussion, however, we will assume that 
(f) v has a narrow width (e.g. the vertical extent of the dots 
defining v s in Fig. IB 1 f) . 

In MCs, and generally in any environment, the velocities 
are not discretized but rather vary continuously throughout 
the medium. Therefore, the velocities in Figure IBTI are more 
accurately described by the dashed line connecting the de- 
fined v s s in each zone. Consequently, at the observed velocity 
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Figure Bl. Velocity component along the LoS, v s vs. s. A sample 
observed frequency (or velocity) is shown by the solid line, u j, s . 
The units of v 3 and s are arbitrary. Due to the large velocity 
gradients, there would not be any observed emission at this i> bsi 
unless the velocities were interpolated (dashed line). 



there would in fact be emission along this LoS, due to the 
overlap of v„ and v b s (at two locations in this example). 

In order to account for the continuous variation in cloud 
properties, we have implemented interp olation into our ra- 
diativ e transfer calculations (also see iPontoppidan et "ail 
2009). All the required quantities, such as temperature, CO 
density, H2 density, and position, are linearly interpolated, 
effectively increasing the physical resolution along the LoS. 

Figure lB2l shows how under-sampling the velocities may 
produce inaccurate line profiles. The figure shows the rele- 
vant quantities from LoS A shown in Figure[l5]and discussed 
in Section f4. 51 The left panels show the results using inter- 
polation (16x the original grid, so As = 0.0049 pc), and the 
right panels show the result using the original simulation 
grid (As = 0.078 pc). Clearly, there are stark differences in 
the spectra. When unresolved, the spectrum has numerous 
"flat tops," and there is a sharp decrease in the intensities 
between —3 and —4 km s _1 . These features do not occur in 
the spectrum obtained from the interpolated LoS. The "fiat 
tops" and intensity drop could be misidentified as satura- 
tion or absorption features, respectively, when in fact they 
are simply features arising due to the under-resolved grid. 

The differences can be understood from the velocity and 
optical depth profiles. In the high density region responsible 
for most of the emission (see Section I4.5J1 , there is a very 
large velocity gradient. At an observed velocity of —4 km 
s _1 , there would not be any emission from this region, due 
to the large gap in the velocities at LoS position ~18.5 pc. 
However, when interpolated (left panels), there are continu- 
ous velocities in this region, resulting in enhanced emission 
at —4 km s _1 . 



For the interpolated 0.0049 pc and original 0.078 pc 
simulations, the integrated intensities from LoS A are 187.1 
and 173.8 K km s _1 , respectively. This results in a slight 
difference of ~10 %, and would lead to a similar discrep- 
ancy in the X factor from this LoS. The precise difference 
would of course depend on the particular LoS — LoSs with 
many regions with large velocity gradients would be strongly 
affected by under-resolved velocities. For model n300, the 
cloud-average X factor is (X)=2.2x 10 20 cm" 2 KT 1 km" 1 
s, regardless of whether interpolation is employed. Thus, for 
the n300 model, only the line shapes along individual LoS 
are affected by under-resolution, but the velocities are suffi- 
ciently resolved for estimating the cloud-averaged integrated 
intensity. 

For the models discussed in Section 4.4 where we have 
directly modified the velocities to have larger dispersions, 
there will of course be large velocity gradients in adjacent 
zones. For these models, we find that interpolation results in 
more significant differences. Recall that in the altered mod- 
els, the velocities are replaced in the original grid utilized in 
the chemo-MHD simulation. The altered velocities on this 
grid have the chosen dispersions. In the radiative transfer 
calculations, due to the large velocity gradients, interpola- 
tion leads to the addition of numerous values to the veloc- 
ity dataset. In the models with large chosen dispersions in 
the original grid, this interpolation can potentially lead to 
a velocity dispersion that is different from the chosen one. 
For instance, in the a = 24 and 240 km s" 1 models, after 
interpolation (with 16 x increased resolution), the effective 
velocity dispersion becomes 21 and 205 km s _1 , respectively. 
In both the a = 24 and 240 km s _1 models, we find that 
interpolation increases the cloud-averaged W by factors of 
2 and 4, respectively, and thereby decreases the X factor by 
a corresponding amount. As these models are highly artifi- 
cial to begin with, we do not ensure that all the velocities 
are well resolved. Nevertheless, qualitatively, we can be sure 
that increasing the velocity dispersion leads to decreasing 
the X factor. As we state in Section 5.3, models of real 
sources with very high dispersions, such as ULIRGs, would 
have to be well resolved in order to quantify the X factor— a 
relationship. 
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